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Aims. We investigate the effect of a geometric inconsistency in the calculation of synthetic spectra of giant stars. 

Methods. Spectra computed with model atmospheres calculated in spherical geometry while using the plane-parallel approximation for line 
formation calculations (s_p), as well as the fully plane-parallel case (p_p), are compared to the consistently spherical case (s_s). 
Results. We present abundance differences for solar metallicity models with T c ff ranging from 4000 to 6500 K and logg from 0.5 to 3.0 [cgs]. 
The effects are smaller for s_p calculations (-0. 1 dex in the worst case) than for the p_p case (up to +0.35 dex for minority species and at most 
-0.04 dex for majority species), both with respect to the s_s case. In the s.p case the differences increase slightly with temperature, while in 
the /?_/? case they show a more complex behaviour. In both cases the effects decrease with increasing log g and increase with equivalent width. 
Conclusions. Within the parameter range of F, G and K giants, consistency seems to be less important than using a spherical model atmosphere. 
The abundance differences due to sphericity effects presented here can be used for error estimation in abundance studies relying on plane-parallel 
modelling. 
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1. Introduction 

Classical abundance analyses are usually based on model at- 
mospheres calculated assuming plane-parallel geometry. This 
means that radiative transfer is solved in only one depth vari- 
able, neglecting the curvature of the atmosphere. For most stars 
this is justified because the extension of the atmosphere is neg- 
ligibly small compared to the stellar radius. Extreme cases are 
cool giants and supergiants, where geometry effects due to the 
thickness of the atmosphere have to be taken into account. They 
have been studied accordingly well in the literature. 

Sphericity effects on structures, low resolution spectra 
and colors of cool (T e g < 4000 K) giants and surfer- 
giants have been studied for example by Plez et 
They showed that for one of their hottest models with 
(r e ff,log^',mass)=(3800K,1.0,l M Q ), the changes in tempera- 
ture structure when releasing the plane-parallel approximation 
are relatively small (about 35 K at logTgp.,., = -4, the atmo- 
spheric extension being about 4%). US (1990) discussed ef- 
fects of atmospheric extension on the determination of effec- 
tive temperature, surface gravity and abundances of cool su- 
pergiants (r e (j < 4500 K). He noted that the effects are not 
linear with temperature and gravity. 

Spherical model atmospheres for a large range of giant 
star temperatures have been presented by lHauschildt et alJ 
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( 1999J). They also presented changes in spectral flux 
with respect to plane-parallel models. For a model with 
(r e ff,log§)=(5600 K,0.0) they show for example that the 
change in central depth of individual spectral lines can be up 
to 15%. 



F to K giants are commonly used to study abundances 
of stellar systems, mainly due to their intrinsic brightness. 
With respect to atmosphere geometry they represent border- 
line cases and are therefore usually analysed using plane- 
parallel radiative transfer. Many examples for recent large scale 
abundance analyses of giant stars can be found in the litera- 
tur e. Cepheid stars at variou s pulsation phases were studied 
bv iLuck & Andrievskvl ll2004h . These stars are used to deter - 
mine the Galactic abundance gradient ( Andrievskv et alJ20 04). 
Anoth er example is the an alysis of probable giant star planet 
hosts ( Taked aet alJ 120051) . Furthermore, red giants are fre- 
quently contained in samples of Galactic halo stars, which are 
used to study Galactic chemical evolution and nucleosynthe- 
sis by means of classical abundance analysis. Earlier examples 
are the large scale high resolu t ion ab undance studies of metal 
poor stars bvlMc William etafl l l 19951 31 stars with logg < 2.5) 
and iRvanetalJ dl996(. 7 stars with logg ^ 2 .0), whose re- 
sults have been used by Travagl io et alJ (J2001) in their study 
of the chemical evolut ion of the Galactic halo. More recently, 
ISimmerer et alJ J2004I) determined abundances of Fe, La and 
Eu for stars over a wide range of metallicities. A substantial 
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fraction of those are giant stars in the logg range of 0.5 to 2.0, 
aU with [Fe/H] < -1 dex. 

In recent years, grids of model atmospheres calculated 
in spherical g eometry have become available to an increas - 
ing extent (e.g. lGustafsson et all2003tlHauschildt et alll9 99). 
However, line formation will probably continue to be computed 
in plane-parallel geometry in abundance analysis procedures. 
In this paper, we investigate the effect of such an inconsis- 
tency on synthetic spectra and abundances of giant stars. We 
then compare to the effects due to using a consistent, but plane- 
parallel combination of model atmosphere and spectrum syn- 
thesis. 

Figure ^ shows the atmospheric extension of spherical 
MARCS models (see Section |3 with £ t = 2 km s _1 , solar 
abundances and various temperatures and gravities as a func- 
tion of bolometric magnitude. Here, the extension is defined as 
the difference in radius at the points in the atmosphere where 
Tross = 10~ 5 and 1 relative to the radius at tr oss = 1. The bolo- 
metric magnitude has been calculated from the model fluxe s 
on the scale defined by IAU Commission 36 (Andersen 1999). 
This figure can be directly compared to Fig. 1 of IPlez et alJ 
ill 992b . It shows that atmospheric extensions increase relatively 
rapidly with effective temperature for logg < 2. One can also 
see small local maxima at about 5 250 K, similar to the ones at 
3200 K found by P lez et alJ dl992>. These features in the shape 
of lines of equal log g can probably be attributed to a decrease 
in partial pressure gradient of diatomic molecules such as CO 
or CN. The same plot for models with abundances correspond- 
ing to halo stars with [Fe/H] = -2 dex looks very similar, ex- 
cept that the aforementioned maxima are not present. This is in 
accordance with the low metallicity where molecules are not as 
important. 

Figure|2]shows how the temperature structures of spherical 
models with solar abundances in the parameter range investi- 
gated here differ from plane-parallel ones. The dilution of the 
radiation field in spherical geometry causes the spherical model 
to be cooler than the plane-parallel one above an optical depth 
of a: 1, and the temperature difference increases outwards. 

The differences seen for the hotter models at large optical 
depths can be attributed to numerical effects in the linearized 
mixing length equations. For these models, the lower boundary 
of the convection zone is very close to the lower boundary of 
the atmosphere. Note that these differences are not visible in the 
spectra and that the differences are less than 1 % (see tempera- 
ture scale for the spherical (r e jf,logg)=(6500 K, 1.5) model at 
the top axis of the figure. 

In the outermost layers, the temperature differences for 
models of equal log g and different T e s show a complex be- 
haviour. In certain cases, models with larger extension show 
less surface cooling with respect to the plane-parallel counter- 
part than less extended models, for example the models with 
log §=1.5 shown in Figure|2 This is very likely due to the effect 
of the change in temperature structure on molecular absorp- 
tion. An inspection of the grid models shows that the changes 
in partial pressures of molecules in the outermost layers, when 
going from plane-parallel to spherical geometry, do not sim- 
ply scale with effective temperature. This applies in particu- 
lar to molecules containing carbon, as well as TiO and SiS. 
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Fig. 1. Atmospheric extension ( r( ^" ^ } - 1) of spherical 
models with £ t = 2 km s and solar abundances as a function 
of bolometric magnitude for M — 1,2, and 5 M e as indicated 
by increasing symbol size. Models with equal logg are joined 
by lines, except where certain T e ff values are missing, and are 
labelled with their \ogg value. Model temperatures range from 
4000 to 7000 K in steps of 250 K for M = 1 M and 1000 K 
for M — 2 and 5 M , increasing from right to left along lines. 
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Fig. 2. Differences between spherical and plane-parallel model 
temperatures (T s - T p ) as a function of Rosseland optical 
depth of models with £ t = 2 km s _1 , solar abundances and 
M - I M Q . Solid (dashed) lines show T eS = 5000 K (6500 K) 
models and line width increases with logg. The scale on 
the top horizontal axis gives the temperature of the spherical 
(r eff ,logg)=(6500 K, 1.5) model in units of 10 3 K. 

Since molecular absorption can eit her cool or heat the outer 
layers ( Gustafsson & Olanderll 19791) . the net effect might well 
be more cooling in a less extended model as compared to a 
more extended model. 

2. Models, spectra and equivalent widths 

We calculated spectra with two different versions of the spec- 
trum synthesis code used in Uppsala (plane-parallel — p — and 
spherical - s) for representative model parameters. We chose 
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Table 1. Number of lines selected for each species in 
the geometrically consistent (s_s) synthetic spectrum with 
(r eff ,logg,£)=(5000 K,0.5,2 km s" 1 ) with equivalent widths 
between 5 and 250 mA, as well as fraction of the element in the 
given ionization stage (in %, averaged over -1.0 < logr« OJJ < 
-0.1). 



Species 


n 


% 


Species 


n 


% 


Species 


n 


% 


All 


9 


0.1 


Gd II 


14 


99.8 


Sc II 


15 


99.9 


CI 


21 


100.0 


La II 


33 


98.6 


Si I 


173 


5.2 


Cal 


28 


0.01 


Mgl 


8 


0.5 


Si II 


4 


94.8 


Cell 


33 


98.3 


Mnl 


17 


0.5 


Sm II 


104 


98.3 


Col 


76 


1.8 


Mn II 


3 


99.5 


Til 


115 


0.1 


CrI 


41 


0.1 


Mo I 


3 


0.3 


Till 


34 


99.9 


CrII 


22 


99.9 


Nd II 


109 


96.7 


VI 


57 


0.1 


Dy II 


5 


99.3 


Nil 


124 


2.6 


VII 


21 


99.9 


ErII 


5 


99.6 


OI 


4 


100.0 


YII 


10 


99.8 


EuII 


7 


98.6 


PrII 


73 


94.8 


Zrl 


6 


0.1 


Fel 


562 


1.0 


SI 


21 


91.1 


ZrH 


7 


99.9 


Fell 


68 


99.0 


Scl 


8 


0.03 









solar metallicity for all models and masses of 1, 2 and 5 M for 
s models. The microturbulence parameter £ t was set to 2 and 
5 km s _I . r e ff was varied between 4000 and 6500 K in steps of 
500 K for 1 M and 1000 K for 2 and 5 M B . The logg values 
ranged from 0.5 to 3.0 in steps of 0.5 for 1 M and were set 
to 1.0, 2.0 and 3.0 for 2 and 5 M B . For 1 M , this corresponds 
to radii at t Ross = 1 of « 5R Q (logg = 3.0), ~ \5R & (logg = 
2.0) and * 90R Q (logg = 0.5); for 2 and 5 M and logg = 2.0, 
the radii are w 20 and 4QR , respectively. For logg = 0.5, only 
models with T e « < 5000 K were available. Normalized syn- 
thetic spectra were calculated for wavelengths b etween 5 400 
and 7 200 A using atomic line data f rom VALD JKunka et a lJ 
1 19991) and MARCS model structures jGustafsson et alJl20 03) 2 
as input. For a discussion of the method for solving the radia- 
tive transfer equati ons in the atmosphe ric mode l s and spectrum 
synthesis see lPlez et all i 19921) and iNordlundl Jl984 . We re- 
gard three different combinations of model atmosphere (atm) 
and spectrum synthesis (syn) geometries: atmsyn = ss (con- 
sistently spherical), s_p (inconsistent), p_p> (consistently plane- 
parallel). 

As a next step, we selected all lines in this spectral re- 
gion which were found to be not blended by more than 30 % 
by neighboring lines. This line selection was done for the 
(5000 K,0.5), £ =2 and 5 km s _1 , 1 M s.s models. Equivalent 
widths were then calculated for these lines for the ss models 
with an equivalent width / abundance fit version of the spec- 
trum synthesis code. Finally, abundance differences were ob- 
tained for i_p and p_p models by fitting the calculated equiva- 
lent widths to the ss ones. Tablef^lists all species contained in 
the line list and the number of lines selected for each of them, 
as well as the fraction of the element in the given ionization 
stage, averaged over -1.0 < logTR OJJ < -0.1. 

2 http://marcs.astro.uu.se 



3. Results 

We take the consistent s_s spectra as a reference and present a 
comparison to the spectra with the two other geometry combi- 
nations. Fig.|5]illustrates the geometry effects on line profiles. 
One can see that mainly the line centers are affected, which can 
be ascribed to the way the model atmosphere structure changes 
(see Fig.|3}. The following two subsections summarize the ef- 
fects on line abundances for 1 M models. We restrict the inves- 
tigation to lines with equivalent widths between 5 and 250 mA. 

3.1. Effects of inconsistent spectrum synthesis 

Figure0]shows the differences between s_p and s_s abundances 
versus s_s equivalent widths (W) for Fe I lines (upper plot) and 
Fe II lines (lower plot). Models with T e g - 5000 K, low mi- 
croturbulence and various gravities are shown. Apart from the 
expected increase towards lower gravities, one notices that the 
effect is on average about 0.01 dex smaller for Fe II than for 
Fe I, for lines with W < 50 mA and log g < 1 . The same dif- 
ference between ionization stages is seen for all other elements 
with at least 10 lines of each of the two ionization stages (Cr, Ti, 
V). The differences are slightly larger for higher r e ffS. This is 
illustrated in Fig. [5] which shows the differences forFe lines for 
£ t = 2 km s~ 1 , log g = 1 .0 and several T e g values, averaged over 
bins of 25 mA. The trends are weaker for £ t = 5 km s _1 and dif- 
ferent for different species. The differences are in general larger 
for the lower £ t value, where line saturation sets in at smaller 
W. Neutral species show the same differences as Fe I, except 
for C I, S I and Si I, for which the effect is on average smaller 
by 0.01 dex for the lowest logg values and weak lines. Singly 
ionized species show the same differences as Fe II, except for 
weak lines of rare earth elements for log g = 0.5. For Ce II, the 
effect is on average smaller by 0.01 dex, and for Gd II, Pr II, 
Nd II, La II, Sm II, the effect is on average larger by 0.01 dex. 

3.2. Effects of consistent plane-parallel spectrum 
synthesis 

Figure [6] shows the differences between pjp and s_s abun- 
dances versus s_s equivalent widths for Fe I and Fe II lines, for 
the models with T e g = 5000 K, low microturbulence and vari- 
ous gravities. In this case, neutral and ionized Fe lines show an 
opposite effect, which is also seen for all other elements. The 
dependence of abundance differences on T e g is more complex 
than in the previous case (Section f3.1> . as shown Figure|4] For 
Fe I lines, the differences are essentially zero at T e ff = 4000 K, 
increase up to T e ff = 5000 K, and decrease slightly towards 
higher r e ff. For Fe II lines, the differences decrease with in- 
creasing r e ff over the whole range. 

Again, most neutral species behave similarly to Fe I, with 
the following exceptions: 

- C I, O I and S I (all predominantly neutral): The effect has 
opposite sign to that of Fe I (similar to Fe II). 

- Si I: The effect is less than half that of Fe I. 

- Ti I and V I: The effect is on average larger than for Fe I. 
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Fig. 3. Comparison of s_s (thick line), s_p (thin line) and p_p (dotted line) spectra for 
(r eff ,logg,£,M)=(5000 K,0.5,2 km s" 1 ,! M ). 



Most singly ionized species show the same differences as 
Fe II. Exceptions are again the rare earth elements, but also Sc 
and Ti: 

- Ce II, La II, Nd II, Sc II: The effect is smaller than for Fe II 
and changes sign for large equivalent widths. 

- Sm II, Ti II: The effect is on average 0.01 dex smaller than 
for Fe II for log g < 1 .0. 

3.3. Mass dependence 

For models with (r eff ,logg,£)=(5000 K, 1.0,2 km s" 1 ) and a 
mass of 2 M , the abundance differences are smaller by 0.005 
(weak lines) to 0.01 dex (strong lines) compared to 1 M mod- 
els. For a model mass of 5 M , they are smaller by 0.01 (weak 
lines and strong lines of majority species in the p_p case) to 
0.03 (strong lines in the s_p case) to 0.05 dex (strong lines of 
minority species in the p_p case). In general, higher masses 
result in smaller effects and the dependence on mass is propor- 
tional to abundance difference. As can be seen in Fig. ^ the 
extensions in 5 M Q , log g- 1.0 models are approximately the 
same as in 1 M , log g=\.5 models. We therefore expect that 
an increase in mass from 1 to 5 M Q has the same effect as in- 
creasing logg by 0.5. This is confirmed by the calculations. 



4. Discussion 

4.1. Trends of effects with line data 

It seems that the temperature structure used as input for the 
spectrum synthesis calculation is more important than the ge- 
ometry used in such a code. Since the temperature differences 
decrease with depth in the atmosphere, spectral lines formed at 
different depths will be affected differently. We calculated aver- 
age formation depths (tr oss ) of the line centers for a few species 
which show different effects (Section [3} for the s model with 
(r eff ,log^,^ t )=(5000 K,0.5,2 km s" 1 ). The distributions (rela- 
tive frequency of lines forming at different log(rR OSS )) are rather 
similar, except for those of C I, S I, Si I and Ti I. They show a 
somewhat higher and narrower peak at log(TR OSS ) * -0.3 than 
that of Fe I. This will introduce a small species to species dif- 
ference in the geometry effects. 

The largest part of the differential effect can however be at- 
tributed to the different distributions of excitation energies over 
lines of the same species (see Section l4.2.3l below). Figures|8] 
and [9] show the correlations between excitation potential and 
abundance difference (for ionized species) as well as between 
the difference of ionization and excitation potentials and abun- 
dance difference (for neutral species). 
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Fig. 4. Abundance differences for inconsistent and consistently spherical equivalent widths of Fe I lines (upper plot) and Fe II 
lines (lower plot) as a function of ss equivalent widths for models with a mass of 1 M Q with £ t =2 km s _1 , 7^=5000 K, and 
different logg values as indicated by different symbols. 



4.2. Analytical estimation of effects 

As giant stars with exactly plane-parallel atmospheres do not 
exist in nature, it will be difficult to verify these results with 
observations. A possible experiment would be to perform an 
abundance analysis of a sample of binary stars, where one 
component is a dwarf star close to the zero-age main se- 
quence and the second components are giant stars of varying 
sizes. When both stars are analysed with a plane-parallel spec- 
trum/atmosphere code, any difference in abundance could be 
attributed to sphericity effects. However, this requires that all 
other physical phenomena, like convection, deviation from LTE 
and scattering, are modelled correctly. Here, in order to conduct 



a "sanity check" on the results, we try to reproduce them using 
an analytical solution to the problem. This will also help to get 
a better understanding for the cause of the effects. Such solu- 
tions require invoking a few more approximations than used in 
the numerical model. We follow IChandrasekharl Jl934l) . who 
derived a solution to the equation of radiative transfer (in polar 
coordinates), including the term containing the dependence of 
intensity I on the angle of inclination to vertical 9 



dl sin dl _ 

cos0- — = tcp(S - I), 

or r o9 



(1) 



where r is the distance from the stellar center (i.e., the ra- 
dius), k is the opacity, p the density and S the source function. 
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Fig. 5. Abundance differences for inconsistent and consistently spherical equivalent widths of Fe I lines (upper plot) and Fe II 
lines (lower plot) as a function of s_s equivalent widths for models with a mass of 1 M Q with £ t =2 km s _1 , log g=1.0, and 
different T e g values as indicated by different symbols. The differences have been averaged over bins of 25 mA, with bin centers 
offset for clarity and error bars showing the standard deviations. 



In local thermodynamic equilibrium (LTE), S is given by the 
Planck function B. Further, if the mean intensity J, the flux nF, 
and the radiation pressure integral K are defined as the zeroth, 
first and second moment of the intensity, respectively, the first 
Eddington approximation 3K = J is used, along with the con- 
dition that at the outer boundary of the star F — 2J (second 
Eddington approximation). Radiative equilibrium is expressed 
by Fr 2 =constant. The dependence of B on optical depth r is 
then 




where ttFr is the flux at r — R, which occurs at t — 0. Note that 
this definition of R is different from the radii given in Section|2 
If tross = 10~ 5 is taken as an approximation for t — 0, R is the 
radius at Tr oss = 1 enlarged by the corresponding extension 
(Fig. ID . We adopt this approach in the following. 

The same approximations as above are subsequently used 
when considering the formation of absorption lines. Now, all 
radiation quantities are indexed with frequency v to indicate 
that they are defined in the range v to v + dv in continuum re- 
gions, and additional primes indicate quantities inside spectral 
lines. 
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Fig. 6. Differences between consistently plane-parallel and spherical abundances versus s_s equivalent widths of Fe I lines (upper 
plot) and Fe II lines (lower plot) for models with a mass of 1 M Q with £ t =2 km s _1 , r e ff=5000 K, and different logg values as 
indicated by different symbols. 



4.2.1 . The ss case 

IChandrasekhaJ lll934h assumed a power law variation of 
with r 



Ty = Cf, 



(3) 



where c and m are constants and m < 0. Furthermore, R was 
assumed to be sufficiently large so that 



1 3 C Ty dr 



R 2 2 



o r- 



2 ' 



Assuming small \m\, he derived a surprisingly simple expres- 
sion for the dependence of normalized line flux on line strength, 



the latter being measured by the ratio of the coefficients of line 
to continuous absorption 77, 

F' 1 



r i 1 (5) 

Fy 1+77 

Equation [5] presumably describes the s_s case and can be 
direc tly compared to a similar expression derived by Eddi ngtor] 
( 1929)' for the plane-parallel case 



F' 



(4) 7T = 



1 + \q 



(6) 



Here, q 2 = 3(1 + 77XI + erf) includes a parameter e measuring 
the importance of collisions (i.e. e = 1 in LTE). Comparing 
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Fig. 7. Differences between consistently plane-parallel and spherical abundances versus s_s equivalent widths of Fe I lines (upper 
plot) and Fe II lines (lower plot) for models with a mass of 1 M Q with £ t =2 km s~', log g=1.0, and different r e ff values as 
indicated by different symbols. The differences have been averaged over bins of 25 mA, with bin centers offset for clarity and 
error bars showing the standard deviations. 



the two expressions, it becomes clear that Eq.[5]by far overes- 
timates the effect of sphericity, since the normalized line flux 
decreases much more rapidly with line strength than in Eq. [6] 
and it approaches zero for large 77, whereas in Eq.[6] the limit- 
ing value is 1/(1 + V3/2) = 0.54. Thus, Eq.|5]cannot be used 
to estimate the effects for the parameter range studied here. An 
inspection of the calculat ed models shows that neither of the 
two assumptions used by ChandrasekhaJ Jl934h are fulfilled. 
We determine values of m and c by fitting the function given 
by Eq.[3]to the calculated values of t 5000 ^ and r in the range 
-5 < log t 50()0 ^ < 0. The ratio of ^ and the right hand side of 



Eq.|4]is then between 0.5 to 0.7 when using r 5000 ^ = 1 as up- 
per integration limit and correspondingly larger for smaller up- 
per integration limits. Furthermore, the values of \m\ are rather 
large (between 100 and 2000). 



4.2.2. The s_p case 



For the s-p case we start from a slightly different expression 
for the normalized line flux, derived for the case where the 
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continuum source function varies linearly with optical depth 
(S v = By = a + bT v ): 

K = b + ag 1 + f V3 

(see lEddingtonll929l his Eq. 22). 

The effect of using a spherical model atmosphere com- 
pared to a plane-parallel one can then be estimated by compar- 
ing corresponding numerical values of the ratio b/a. In LTE, 
S v - By - C v /(e x - 1), where Cy is constant for a given fre- 
quency and x = hv/kT, T being the temperature and /; and k 
the Planck and Boltzmann constants, respectively. The deriva- 



tive of the logarithm of the Source function with respect to t v 
is then 

1 dS v xe x 1 dT 
S v dr v e x - IT dr v 

For a grey atmosphere, where the continuous absorption is in- 
dependent of v, the temperature structure is given by T 4 cc 
(1 + It) in the plane-parallel case. Inserting in Eq.|8] together 
with the linear source function law, we derive b/a = |x at 
t = 0. This expression only holds for x > 3, for which 
e x /(e x - 1) can be approximated by 1. Note that the lowest 
value of x occurring in the atmospheric models studied here 
is about 3 (for 7200 A and the hottest models at tr oss = 1). In 
the spherical case, using Eq. |2j the grey temperature structure 



10 



U. Heiter and K. Eriksson: Geometry test 



x 

CD 
"D 



CD 
CJ 
C 
CD 



CD 
O 
£= 

cC 
"D 
C 
=3 

< 



0.40 
0.35 



-g 0.30 



§ 0.25 

CD 
s 

CD 

£ 0.20 

T3 
CD 

I 0.15 

CO 
"D 
CZ 

E o.io 
< 



0.05 
0.00 



pp - ss, logg = 0.5, ^ = 2, T 



eff : 



5000 



Crl 
Ti I 
Fe I 
Ni I 
Si I 



y « .1 




3 4 5 6 

Ionization - excitation potential [eV] 



0.04 



0.02 



0.00 



-0.02 



-0.04 



pp - ss, logg = 0.5, ^ = 2, T 



eff : 



5000 




logx R =-0.5 



logx R =-0, 
logg = 0.5 



-0.06 

1 2 3 4 5 6 7 

Excitation potential % [eV] 

Fig. 9. Abundance differences p_p - s_s versus difference of ionization and excitation potentials (neutral species, upper plot) and 
versus excitation potential (singly ionized species, lower plot). The lines in the lower plot indicate the variation of abundance 
difference as estimated from the difference in model temperature only (see text, Section l4.2.3> at log tr OS s = -0.5 (dashed) and 
log tr oss = -0.1 (solid), for different log g values as indicated by different line width. 



can be written as T 4 cc (l + | yrdTj. We use the same law 
for r(r) as above (Eq. |3j to evaluate the integral in the temper- 
ature structure equation and further the derivative in Eq.|8] For 
large values of \m\, this leads to bja = ^xc 2 ^ m R 2 . In order to 
estimate numerical values of the discriminating factor c 2 ^'"R 2 , 
we resort to determining them from the atmospheric models: 
for logg values of 0.5, 1.5 and 2.5, typical values of the factor 
are 1.2, 1.1, and 1.02, respectively. The resulting variation in 
normalized line flux is shown in Fig. ^| as a function of line 
strength rj, for a value of x = 4 (corresponding to e.g. a wave- 
length of 6000 A at a temperature of 6000 K). Also shown is 



the relation for the p_p case (Eq. [6] which derives from Eq. 
when x = 4). Comparing this figure to Fig.[5]shows that the an- 
alytic estimation of the sphericity effect agrees fairly well with 
the modelling outcome in this case, at least for lines of species 
in the minority ionization stage. 

The results of this discussion suggest another observational 
test. When central depths of saturated lines are measured in gi- 
ant stars of different radii, they are relatively unaffected by star- 
to-star abundance differences. Those lines having cores that 
form close to LTE should then be deeper for the more extended 
stars than for the less extended stars by the amount predicted. 
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Fig. 10. Normalized flux as a function of line strength show- 
ing the s_p effect as estimated according to Eq. where 
b/a = §jec 2/m /? 2 , for x = 4 and three different factors c 2/m R 2 
determined from the atmospheric models for different log g val- 
ues. 



4.2.3. Dependence on excitation potential 

The variation of the abundance differences with excitation po- 
tential (x) can under certain circumstances be directly esti- 
mated from model temperature differences. We use the assump- 
tion that for a weak line of a given species, changes in logarith- 
mic abundance are equivalent to cha nges in 5040 /T -y for a line 
from the majority species (see e.g. lGravll 19921 Chapter 14). 
The temperature difference in the region of log tr oss between 
-0.1 and -0.5, where most of the lines form, is about 30, 15 
and 5 K for models with logg = 0.5, 1.5 and 2.5, respectively. 
The corresponding variations of abundance difference with x 
are shown in Fig.|9](lower plot) as lines and are comparable to 
the outcome of the detailed models. For lines of species in the 
minority ionization stage (all neutral elements except C I, O I 
and S I, see Table 0, the differences will also depend on the 
ionization potential. 

4.3. Implications for abundance analysis 

The systematically different effects on lines of different ion- 
ization stages of the same element imply an effect on ioniza- 
tion equilibrium, which is sometimes used to infer gravities 
of pr ogram stars. Just as an exam ple, we take the star YZ Sgr 
from lLuck & Andrievskvl J2004 . At phase 0.249 it has a tem- 
perature of 5500 K, logg = 1.4 and £ t = 3.7 km s _1 . The 
av erage abundance difference (p-p - s_s) for Fe I lines used 
bv lLuck & Andrievskvl (l2004h . which are also in our line list, 
amounts to +0.03 ± 0.01 dex (for the 5500 K/1.5 model with £ 
= 2 or5 km s _1 ). ForFe II, the difference is -0.010±0.003 dex. 
The same difference is obtained by decreasing log g in the s_s 
models by about 0. 1 dex or increasing the temperature by about 
40 K. These valu es are equal to the f ormal uncertainty given 
for this star. In the Tak eda et alJ J2005I) sample, stars at 5000 K 
have logg > 2.5. If we take as an example the star HD 192944 
(£t = 1.4 km s _1 ), the abundance differences (p_p-s_s) are only 



+0.009 ± 0.006 dex and -0.008 + 0.004 dex for Fe I and Fe II, 
respectively. They can safely be ignored until abundance anal- 
yses reach significantly higher internal accuracies compared to 
today. Abundance differences for individual lines and/or stellar 
parameters can be obtained on request from the authors. 

5. Conclusions 

1 . We recommend the use of spherical model atmospheres in 
abundance analyses for logg < 2 and 4000 K < T e ff < 
6500 K. Thus one can avoid both systematic errors on abun- 
dances and differential effects, which can lead to additional 
uncertainties in stellar parameters. Alternatively, restricting 
the analysis to weak lines (W < 100 mA) and/or major- 
ity species will minimize the uncertainties on abundances 
caused by assuming a plane-parallel geometry. 

2. Geometry has a smaller effect on line formation than on 
model atmosphere structure. It is more important to use 
a spherical model atmosphere than to be consistent when 
calculating the spectral lines, i.e. an s_p calculation is bet- 
ter than p_p if tools for s_s calculations are not avail- 
able. Note, however, that this holds only for the range 
of parameters studied here. For a model of a cool gi- 
ant wi th (!T e ff,logg)=(2800 K,-0.6) for example, Aringer 
J2005I) showed that an inconsistent treatment of geometry 
(s-p) results in larger differences in spectral line depth than 
for a fully plane-parallel calculation. 

3. For both cases (s_p and p-p), the results presented in 
Section [3] are intended to provide a guide to the estima- 
tion of systematic errors in abundance analyses due to ge- 
ometry effects. For s_p calculations, the largest differences 
encountered with respect to the s_s case are -0.1 dex. In 
the p_p case, differences are up to +0.35 dex for minor- 
ity species and at most -0.04 dex for majority species. The 
abundance studies of halo stars mentioned in Section^con- 
tain stars in the log g range where sphericity effects are im- 
portant. Although the abundance difference values given in 
Section|3]are based on solar abundance models, effects are 
expected to be similar in low-metallicity models, since the 
atmospheric extensions are similar. 
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